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Abstract 

II ' In this paper hadroproduction of charmonium states in polarized pp collisions is dis- 

Mh. cussed. A thermal picture for the gluonic cloud of comovers is given making contact 

r~| ! between the formalism and the measured unpolarized cross sections. The experimentally 

observed non-polarization of the final J/'i/; states leads to the consequence that no corre- 

. lations between the initial proton spin and the final charmonium spin should be existent. 

r> I Hence the single spin asymmetries vanish to leading order in that model. 



PACS: 12.39.Hg, 29.27.Hj, 29.25.Pj 



1 Introduction 

Charmonium production in polarized pp scattering at RHIC is one of the key experiments to pin 
down the polarized gluon distribution amplitude [|l|. The standard formalism in which heavy 
quarkonium production in hadron hadron collisions is described is still the Color Octet Mecha- 
nism (COM) [0]. However, it is known that this mechanism fails to describe the experimentally 
observed non-polarization of the final J/ip and ip' meson [Q and, furthermore, it predicts a 
Xci/Xc2 production ratio which is far too low 0]. Recently, it has been shown that these prob- 
lems can be cured by assuming that the charmonium formation happens through rescattering 
with a gluon cloud of hard comovers p, ^ : The two colliding hadrons form through the self in- 
teracting gluon field a gluonic medium in which the heavy quarkonium formation is directed by 
hard rescattering processes. This means a crucial qualitative difference as to electroproduction 
where one of the collision partners is a lepton and cannot participate in strong interactions. 
As a deeper understanding of the heavy quarkonium formation mechanism in polarized pp 
scattering will be very important for the extraction of the polarized gluon density, we want to 
investigate what consequences this theory implies for the charmonium production in polarized 
pp scattering. In addition to the rescattering picture developed in ||, ^ we will present a ther- 
mal description of the comover cloud. 

Double spin asymmetries in charmonium hadroproduction have been calculated in the frame- 
work of the COM formalism, which is based on a systematic non-relativistic velocity expansion 
101, for the prospective HERA-iV experiment 0, ^ |10|, |TT|, [1^. In the framework of the Color 



Singlet Model (CSM) [|^, |T4| double spin asymmetries in J/i/: production have been studied 
in [|l^ . The CSM has been shown to be incompatible with the absolute size of the unpolarized 
J/ip cross section [^]. Attempts to cure this failure by the fc^-factorization approach within 



the CSM have not led to satisfying results [16, O]. Therefore one is looking for a suitable 



combination of the CSM and the COM formalism in the fc^-factorization approach which may 
release the polarization problem of the final J/ip and ip' [|18| . A discussion for the double spin 



asymmetry for RHIC energies in terms of the CSM formalism can be found in |]19| . Higher order 
velocity corrections in the polarized case in terms of the COM formalism have been discussed 
in [^, |21], ^. It is the aim of this paper to add to this discussion the possible insights the 
comover rescattering picture can offer as to polarized hadroproduction of charmonium states. 

In Sees. ^ and ^ we will write down the total cross section for polarized S-wave and P-wave 
charmonium production. For the gluonic comovers we will develop a thermal description as a 
boson gas and explain the consequences for the polarized partonic cross sections, in Sees. § 
and 1^ we will fit the parameters of the theory (volume of the comover cloud, the energy the 
charm quarks carry on the average, and the expansion parameter p of the charmonium system) 
to the numerical values of the measured unpolarized charmonium cross sections. We will dis- 
cuss the implications on the energy transfer from the colliding particles to the cloud and from 
the comover cloud to the charmonium system. We give also a picture of the geometry of the 
cloud. Finally, in Sec. ^ we will calculate the double spin asymmetries for inclusive J/ip, ip' 
and XcJ, J = 0,1,2 production in the framework of the thermal description. 




J/y, y' 



Figure 1: Typical diagram for the S-wave charmonium production amplitude. The formation 
of the J /ijj mesons happens through hard interaction with a gluonic cloud of comovers (F). 



2 Cross section for S-wave quarkonium production 

The gg- fusion amplitude for ^Si quarkonium production in the presence of a background gluon 
field F^(^) (see Fig. |^ for a typical diagram) has been derived in [^ Q: 
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Here Rq is the value of the S-wave quarkonium wave function at the origin. 2m is the energy 
of the two incoming gluons. To the accuracy of the first order velocity expansion used here it is 
identical to the energy of the charm quarks before and after the interaction with the comover 
cloud and also identical to the mass of the J /'ip or ip' produced. Effectively, we will see later 
that we have to adjust the numerical value for m to experimental data of the unpolarized 
cross section. Ai, A2 are the helicities of the two incoming partonic gluons with polarization 
e(Aj), z = 1, 2. It has been pointed out in [^ ^ that the fact that the polarization of J /if) and ip' 
is small and consistent with zero, as can be seen from fixed target data |2^, |2^, |25|, ^, ^ ^ 
and also from the new CDF data 0, leads to the condition that in the cm. system of the two 
quarks the relation: 

|Fg(£)|2<|r^(£)x£|V^' (2) 

must be fulfilled. F^ = —j==V^ is the four-vector potential of the interacting gluonic field. 

The philosophy of our concept is that the quarkonium production happens in a heat bath of 
essentially co-moving (real) gluons. The polarization tensor of the real gluons can then be 
described by: 

r,(^)r;,(£) = -g,, + ^^^^'/^^^'^^ - ^#^ . (3) 
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The introduction of a real gluon field means on the other hand that the incoming charm quark 
from the gluon gluon fusion amplitude must be slightly off-shell. In the calculation this off- 
shellness is not included. In fact, it can be regarded as a velocity correction of order v = 2lQ/m 
(see App. ^ which we will neglect in the calculations that follow. The other alternative, i.e. to 



have a virtual field F leads together with the requirement Fq ~ to unphysical consequences 
as we show in App. |B|, whereas all requirements turn out to be natural for a real gluon field F. 
For the ip' and Xci asymmetries discussed in Sec. |] all model dependent parameters cancel and 
we coincide exactly with the model for F proposed in p, Q. 

In order to fulfill the condition that ensures the non-polarization of quarkonium states Eq. (j^), 
we parametrize n = (1, be), with e being a unit vector and < 6 -C 1. In fact, for 6 = we find 
Fq = 0. One can interpret b approximately as the relative velocity of the gluon comovers to 
the cm. system of the two incoming gluons and e as the direction of this relative movement as 
to the axis defined by the two incoming gluons which is the z-axis in our case. In the following 
we will refer to b then as displacement parameter. 

The heat bath, in which the quarkonium production takes place has the temperature T = l/j3 
and is described by the Bose-Einstein statistics. This takes into account the multiple gluon 
interaction inside the cloud. Guided by this philosophy, we define the following three quantities: 
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One should notice that with this form we reproduce the Stefan-Boltzmann law, that the energy 
density represented by the field squared grows ~ T^ with temperature. The polarization P can 
take the values P = —1, 0, +1, and the polarization vector e(P), (and also g{Sz)) is defined by: 



r (-i,-2,o)/v^ if p = +i 

e(P) = <^ ( 0, 0, 1) if P = 

i (+l,-z,0)/v^ if P = -l 

We then can write down the following partonic cross section: 



(5) 



27r 



C^AiA2S^ 



E 



2(2m)3 ^ (iV2 - 1) 



5 Ti^alRl 



M{^Si,S, 



9 (2m) 



-Xi67^' (<5^ Ff *i - <5^^F§°1 - <5^^f5°*1, 

J^Ai \ Jz Jz-^yj Jz 'JzT^2 >^z Jzi^l 



5 -n^alRl 



9 (2m)6 

The corresponding hadronic cross section is given by: 
^ dcrhAiA2S:, 



5-x.'' (ri? + rLIJ + K^'stri^.' + Ai^r.'^ (i - 25l) r^ 
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Now one can isolate the various components: 
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where (2m,) is the energy of the two incoming gluons. Eq. (|^) is vahd in general, even without 
any assumptions about the thermal nature of the gluonic cloud, if we modify the definition 



of r*^^'-^'™*^ accordingly. In case we apply our model we find always r^°_^^ = 0. Then, the 

single spin asymmetries in our model are proportional to SzT^J , which means proportional 6^. 
As we know from the non-polarization of the final J/ip and ip' meson that b must be small 
and consistent with zero, we consequently predict the absence of single spin asymmetries in 
charmonium hadropro duct ion. Furthermore, if we can set b to zero as the non-polarization of 
the final J/ip and ip' suggests, we see that no initial spin - final spin correlations are present 
and the double spin asymmetry reduces to: 

Soo(5.) ^ G(xi,(2m)2)G(x2,(2m)2) " ^ ^ 

In other words, the experimental fact that the J/ip and ip' are produced in a non-polarized 
mode assures that the double spin asymmetry Eq. ( |T0| ) is essentially only dependent on the 
ratio of the unpolarized and polarized gluon parton distributions, which is a very important 
statement as far as the possibility is concerned to isolate the polarized gluon density from direct 
J/ip or ip' production data. 

We can describe the gluon cloud of comovers with a Bose distribution of momenta in the 
sense of a gluon plasma, in which the formation of the charmonium states takes place. To get 
a quantitative model the following quantities have to be set: 

• The temperature T of the cloud which should lie above the phase transition of hadronic 
matter, i.e. larger than typically 200 MeV. It should also be much smaller than the 
typical charm mass of 1.5 GeV, otherwise the interaction with the gluon cloud would 
rather inhibit the charmonium production than catalyze it. 

• The volume V of the cloud which should be large enough to comprise the cc system. 

• The energy (2m) of the two initial gluons. 

• The expansion parameter p from the quarkonium wave function at the origin to the 
interaction point with the comovers. 

From the non-polarization of the final J/ip and ip' we can already set the parameter 6 = for 
the following. We will do this analysis in Sec. || and Sec. |^ after having collected the cross 
sections for the P-wave charmonium production in the next section. 

3 Cross section for P-wave charmonium production 

P-wave charmonium production for the mesons Xco and Xc2 can occur via two gluons without 
any contribution from comovers, see Fig. 0(a). This corresponds then to the contribution being 
calculated from the Color Singlet Model (CSM). The comover contribution for the production 
of Xcj, J = 0,1,2 mesons becomes important only at the C(«g) level, see Fig. 0(b,c), where the 
comover contribution is supposed to be dominant over all other CSM contribution of the same 
order in a^. The 0{a'j.) CSM partonic cross sections are given by: 
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Figure 2: Typical diagrams for the P-wave charmonium production: (a) CSM contribution 
without comovers at (9(a^); contributions with comover interaction at 0{al): gg fusion (b) 
and qq annihilation (c). 
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We show the derivation of these cross sections in App. ^ to make also sure for the constants 
needed to make contact between the thermal amplitudes and the physical cross sections. Av- 
eraging over the gluon helicities Ai and A2 and summing over all possible final states Jz we 
reproduce the unpolarized partonic cross sections as given in |^9|, ^, ^. At the level of three 
gluons, i.e. 0{al), we can again use the formulas derived from hard comover scattering. This 
time, however, it is hot possible to include a finite displacement b as it would require a further 
expansion of the gluon fusion amplitude, taking into account the Lorentz transformation from 
the gluon gluon cm. to the final quark quark cm. 0. Such an expansion would destroy the 
simple structure of the theory developed so far and therefore we leave it out here. Furthermore, 
as data indicate from the non-polarization of the final J/ip and ip' state should b be small and 
consistent with zero. With this assumption the P-wave amplitude reads O: 
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We define now additionally: 



■■- j=o,o — ■'-0 



^(V) ^ f dH 3 |[e*(J,) X (r(£) X £)J-r 3 C(4) / p 



dH 3 |[e*(J,)x(r(f)xf)rr 3C(4) 

2|£|(27r)3 2 exp(/3(n ■ £)) - 1 2 7^^(3^^^ ""^^ 

.(V) ^ ,./ dH Je3.(^.)rox£rf _ 3C(4) /4,o , ,^, 
^=''-^^ 7 2|£|(27r)3 exp(/3(n ■ £)) - 1 2 7r2/3n3 



^of^ o^l+^tM • (13) 



The explicit form of the tensor e%^^{Jz) necessary for the calculation performed here can be 
found in [^]. We find then for the partonic cross section of the contribution from the diagrams 
of type Fig. g(b): 

(2m) 

This formula states that there is no correlation between the helicities Ai, A2 and the total angular 
momentum J^Jz- P wave quarkonium production can also occur through qq annihilation, see 
Fig. 0(c). Using the formulas given in [^], the amplitude for this case is given by: 
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e*(J,)xe(Ai + A2))(2mrg£ + £2r-) ( J = 1) . (15) 
(^ = 2) 



Here Ai, A2 are the helicities of the incoming quark and antiquark, respectively, i, j are their 
color indices and a is the color index of the gluon from the heat bath. Then the partonic cross 
section reads: 

= 5H^f^r!S)(A, + A,). (16) 

3(2m)io ^^.v 1 2; K J 

If we refer to our model, discussed in the previous section and neglect Fq we obtain: 

rfi(Ai + A2) = (l + |Ai + A2|)r^^ |(l-|e(Ai + A2)-e(J.)p) ( J = 1) .(17) 

"" ^ [ (J = 2) 

The essential statement is that the contribution (c) contains a correlation between the spin of 
the charmonium and the spin of the two incoming gluons. On the other hand this contribution 
(c) is suppressed by a factor T^/(2m)^ versus the contribution (6). As T^/(2m)^ ~ |£p/(2m)^ 
we can again neglect this contribution to the accuracy of the calculation as it is of the same 
order as the velocity corrections not taken into account here. 

The formalism used to derive the partonic cross sections is NRQCD to leading order. In 
principle to the order of accuracy of this approximation the masses of all charmonium mesons 



are the same , i.e. Mjj^ = M^i = M^^^ j=o,i,2 = 2mc 0]. In the same order of accuracy we can 
also identify m = rric in all the formulas above. One has to keep in mind that in such a crude 
approximation the contribution of velocity corrections may be quite substantial. Unfortunately, 
the inclusion of velocity corrections will destroy the simple structure of the relations derived 
in p, 1^ and is therefore out of the scope of the discussion here. To face this problem we will 
fit as a pragmatic ansatz m for the different charmonium states involved in a way that most 
experimental facts are reproduced, m plays then the role of the average gluon energy involved 
in the production of the charmonium state under consideration. 

4 The choice of the quark-energy m 

The starting point of the consideration is the direct J/tp cross section. Here we can as a basis 
identify 2m = Mj/^, because the difference between the Mj/^ = 3.097 GeV and 2171^ = 3 GeV 
is small: 



= J dx,dX2Y.J:ifJ^\S,) . 



3 + ^(ll + 2e^) 



To the order of accuracy of expansion we could apply the same arguments also to the direct ip' 
production. However, such an approximation, which would be in accordance of the standard 
COM velocity expansion, where all charmonium masses are identical Mj/^= ML = M^cj = 
2mc, is in practice very unsuitable as the masses enter in high powers in the partonic cross 
section. We will therefore use an effective mass for the direct ip' production: 



eff^^2^4 



^dir.(V^') = J J dxidx2J's{Ml,,s)G{xuMl,^s)G{x2,M^, 
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.(19) 



We can then fit M^/gg ||32| to the measured ratio <j{i'')dh./o'{J/'ip)dh. given in |^. It has to be 



noticed that this ratio is completely independent of the cloud parameters. The result of the fit 



is: 



M^/eff = 3.4228 GeV , (20) 

which is a bit smaller than the real ip' mass of M^/ = 3.6860 GeV. It indicates in our language 
that there is a net transfer of energy from the gluon cloud into the quarkonium system through 
hard comover rescattering or to stay in a thermal picture that there is a transfer of energy from 
the hotter gluon cloud to the colder charmonium system. The result of the fit can be seen in 
Tab. |1|. 

We can now pay attention to the Xcj niasses. As no comovers enter in the GSM contri- 
bution we should use in this case the original Xcj masses, i.e.: 

^^^'''''Hxcj) = \ E Jdxidx2G{xuMlJG{x2,Miy,^\fPj,J^)6(l-^] . (21) 
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For the contribution resulting from comover rescattering we have to fit two parameters, first 
the effective XcJ niasses and then also the expansion parameter p. We will proceed as follows. 
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0-(V'')dir./0"(-^/^)dir. 


GRV 


CTEQ5L 


experiment E705 


pA 

ttA 


0.23 
0.23 


0.20 
0.21 


0.21 ± 0.05 
0.23 ± 0.05 



Table 1: Fit results using the best value M^i^s = 3.4228 GeV in comparison with experimental 
data from E705 ||32|. The pA experiment was done aX E = 300 GeV and the ttA experiment at 
E = 185 GeV. For the fit we get x^ = 0-08 P^r degree of freedom. 



We take for all three Xcj the same effective Xcj rnass in the spirit that the effective mass should 
be lowered proportional to what was the case for the ip' particle: 
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M,// 



3.2451 GeV . 



(22) 



For the contribution resulting from comover rescattering we can make a fit of the expansion 
parameter p to the E705 data by considering the reduced x fraction, which is defined in a way 
that it is independent of the gluon cloud parameters in our approach: 
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(23) 



In the framework of our thermal description this ratio is independent of the gluon cloud pa- 
rameters. Unfortunately, the reduced x fraction is not directly measured which brings in an 
extra dependence on the parton distributions used. The fit to the reduced ^-fraction yields: 



p = 4.40 



(24) 



This value is a bit larger than the value used in Ref. 0, where the effective mass M^^jcs was 
set to be 2mc- This points to the problem that the relative big expansion parameter means 
that it may be inconsistent to consider only the wave function of the quarkonium system at 
the origin. Besides the substantial velocity corrections this is the second indication that the 
NRQCD approach in general is not a suitable description of the problem. In fact, future analysis 
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(X - frac)^^d 


pA 
experiment 


fit 


ttA 
experiment 


fit 


GRV 
CTEQ5L 


0.29 ± 0.04 
0.30 ± 0.04 


0.34 
0.33 


0.36 ± 0.03 
0.36 ± 0.03 


0.34 
0.34 



Table 2: Fit results using the best value p = 4.40 with Myi ^s = 3.2451 GeV in comparison 



with the values extracted from E705 experiment |^ 

300 GeV and the irA experiment at E 

freedom. 



The pA experiment was done at £" = 
185 GeV. For the fit we get x^ = 0.75 per degree of 



will have to find ways to go beyond the NRQCD approach which we have followed here to be 
compatible with the standard literature of the field. For the details of the fit using p = 4.40 
and M^^^eff = 3.2451 GeV see Tab. |. 



5 Determination of the other parameters of the theory 

The two parameters that are left undetermined so far are the active Volume V of the gluon 
cloud and its temperature T. From the design of the theory the temperature is limited within 
tight bounds. It must be well above Aqcd in order to make the interaction hard and to justify 
the use of perturbation theory. On the other hand it must be smaller than the mass of the charm 
quark rric, otherwise the comover interaction would rather destroy the charmonium system than 
to catalyze it. Now we can make an ansatz taking a constant value T = 500 MeV and let us 
check now what consequences this has for the active volume V. For this purpose we will fit V 
to the data available for inclusive J/ip and ip' production. The situation is simple in the case 
of the ip' production because it is a purely direct process: 
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(25) 

The case of J /ip production is more complicated because to a considerable amount J/ip mesons 
can be produced indirectly via the decay of XcJi </ = 0, 1, 2 mesons predominantly through pho- 
ton emission, [xcj -^ J /"^^ + 7) • Therefore we have to write for the total inclusive unpolarized 
cross section for J/ip hadroproduction: 

cTincxiJ/^) = cTdirXJ/^) + E ^(Xcj)Br(xcj -> J/^ + X) + ddir. (^')Br(V^' ^ J/^P + X) . (26) 

J 

The total Xcj cross section is then given by the contribution from the diagrams in Fig. Q (a) 
and (b): 
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Numerically, the following branching ratios are used [33 



Br(xco^W + 7) = ( 6.6 ±1.8) X 10^3 

Br(xci^W + 7) = (27.3 ±1.6)% 

Br(xc2^W + 7) = (13.5 ±1.1)% 

Bt{^' -^ J/^ + X) = (54.2 ±3.0)%. 

So putting all components together we get for the total inclusive J/ip cross section: 
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For the data of the total cross section n + A 



basically to ||3^ and add the more recent values from [^ 



J/ip, ip' and p(p) + A 
61^, 



(30) 



J/ip^ip' we refer 



27 . The cross sections 



have been rescaled to give the value over the whole range of xp (The details as to this rescaling 
are explained later in this chapter). We reproduce essentially the figures in |Q, except for the 
data point from ||2^ for the nN —>■ ip' cross section, which is displayed a factor 2 too small. 
For the gluon parton distribution of the proton we use the two leading order sets from CTEQ5 
P^ and GRV98 [^]. For the gluon parton distribution in the pion we use the leading order 
parameterization given in |^ (GRS99). For a^ we use the one-loop formula: 



as{f/) 



An 



,(4) 
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Uf) In (^V(ASi) 



using A^^D = 200 MeV , 



(31) 



which comes close to the value used in the GRV and CTEQ5L (leading order) gluon distribution. 
The scale /x^ is given by the only scale relevant for the partonic subprocess of quarkonium 
production, i.e. the quarkonium mass, so /i^ = M'^ 
at the origin Rq is determined to leading order by the decay to e"*"e 



M|,, etc. The quarkonium wave function 



r(j/^,^' 



e e 



4^elaLRl 



(32) 
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Figure 3: The active volume of the gluon plasma fitted from J/^p and ip' data in pA collisions 
(right) and ttA collisions (left). 



We take the values T{J/ip 



MeV, M 



?/i' 



3686 MeV, 



ftp 



number. In order to fix \RW 

data basis here is not very conclusive p3 

and 



e+e") = 5.2374 keV, T{ip' -^ e+e") = 2.3545 keV, Mj/^ = 3097 
a = 1/137 |3^. Cc = 2/3 is the charm quark charge quantum 
we could try to extract it from the decay Xcj -^ 11, however the 
and, therefore, we take here in accordance with [^ 

i.e.: 



the value resulting from the Buchmiiller-Tye potential given in ET 



I p' I 



0.075 GeV^ 



(33) 



It should be noticed that other potentials, also tabulated in [^] yield considerable larger values 
for l-R'iP, up to nearly a factor of 2. This means that the expansion parameter then will be 
reduced by a factor v2 which will not help as to the principle problem mentioned above. 



Using Eqs. (pSf ) and (|30|) the active volume V can be fitted to the data available from pA 
and irA collisions. The result is shown in Fig. ^. We assume a linear dependence in the double 
logarithmic scale, i.e., a dependency of the form V = c f s/GeV^j . The numerical results of 
the fit are shown in Tab. ^ 

It is now the place to make a few statements as to the physical meaning of V and its ge- 
ometry. V is the size of the gluon cloud at the moment the interaction with the charmonium 
pair takes place. It should decrease with s as the faster the collision happens the less time 
has the cloud to form and to expand. Fig. ^ shows the geometry of the cloud. The majority 
of all J/ip are produced at small p±, so we can think us a situation where the cc pair moves 
essentially in the beam axis. In transversal direction the cloud should have a radius roughly 
comparable to I/Iq = l/\i\ ~ 1/7", its longitudinal length however depends on a momentum 
Pll = tt/{VT'^). For simplicity we did not take into account this geometry at the 'thermal' 
integration d^i in Eq. (H) etc. p\\ will grow with s. The ratio Xdoud = P||/v^ which is displayed 
in Fig. ^, shows what fraction of the energy of the system is transferred to the cloud. Whereas 
the CTEQ5L gluon distribution leads to a rising fraction (which is rather unphysical), the GRV 
gluon distribution predicts more or less a constant fraction Xcioud ~ 0.5% for pA collisions and 
a^cioud ~ 1% for ttA collisions. 
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c [fm^] 


GRV 

V 


x' 


c [fm^] 


CTEQ5L 

V 


x' 


pA 

■kA 


5.9212 xlO° 
2.6175 xlOi 


-3.7885x10-1 
-5.9388x10-1 


6.5601 
3.8535 


2.1032 xlO^ 
1.6390 xl02 


-8.2425x10-1 
-8.5073x10-1 


5.7948 
4.2729 



Table 3: Numerical results for the fit of the active volume of the gluon cloud displayed in Fig. ^ 
For the fit a functional form V = c ( s/GeV^j is assumed. 
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Figure 4: Geometry of the cloud of hard rescattering comovers. 



We are now in a position to predict what cross section we will get with our fit at higher 
energies and especially with RHIC energies with s/s = 200 GeV. Fig. ^ shows the fit we made 
in terms of the total inclusive J/ip and ip' cross section. The data have been rescaled to use 
the full range in x^? G [—1,1]. For the tcA collisions we have assumed an xp distribution of the 
form da/dxp ~ (1 — |0.18 — xf\Y, with c = 2.5 for J/ip and c = 3.9 for tp' [^. In case of the 
pA collisions a symmetric xp distribution is assumed. It turns out that for the latter one the 
CTEQ5L distribution predicts an unphysical decreasing cross section at large ^/s therefore we 
will not consider this distribution in the following any longer, whereas GRV shows in all cases 
a reasonable relaxing rising behavior. We can now have a more closer look at the details of the 
various processes contributing using the GRV set only. The first important quantity of interest 
in the Xi/X2 — ratio. It is defined by: 



Xl/X2 



ratio 






(34) 



A comparison of the data and our results is shown in Tab. |. It turns out that the values for 
the Xi/X2 ~ ratio are a bit smaller, but well inside the error bars of the measured irA reactions. 
They would be reproduced even better, if we neglected the GSM contributions altogether. In 
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Figure 5: The fraction of energy XdonA transferred from the projectiles to the gluon cloud in 
pA collisions (right) and it A collisions (left). 



Reference 


E beam [GeV] 


(Xl/X2)exp. 


(Xl/X2)thcor. 


E705 p2| 
E506 35 


185 
515 


1.4 ± 0.4 
1.2 ± 0.4 


1.0456 
0.9194 



Table 4: Comparison of the measured Xxjxi ~ ratio in -nA reaction with the result of the 
theory described here. For the gluon distribution function the GRV98 x GRS99 set is used. 

this case we get (xi/X2)theor. = 1-2 independent of the beam energy, if we set all masses equal 
to Ivdc- So, in principle the theory is capable to describe the large X1/X2 — ratio found experi- 
mentally in contrast to the standard COM and CSM theory. 

Fig. ^ shows some details as to the various subprocesses contributing to the inclusive J/?/' 
production in "pA and vrA scattering. It is shown that the fraction of directly produced J jil) 
versus the whole inclusive cross section decreases with increasing energy until it falls down to 
about 40% at RHIC energies. The s dependence of the cr{ilj')/a{J /ip)^^,: is governed by the scale 
dependence of the gluon parton distribution (here GRV) involved. It goes towards a constant 
for large s which is about 1/4. The XcilXci ratio is rapidly falling with s. In all cases the pA 
and 11 A curves behave similar. 

6 Double spin asymmetries 

The discussion on the partonic cross section has shown two important consequences: 

• As long as we can set the displacement parameter h to zero, which in accordance to the 
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Figure 6: Total cross section for S-wave charmonium production a in nb versus ^/s■. On the left 
the cross sections for J /ip production are displayed and on the right the ones for ip' production. 
The upper two figures show the results for pA and the lower two figures for t^A collisions. All 
nuclear effects have been rescaled so that in principle all the cross sections should display the 
result for pp and 7rp collisions, respectively. The solid and the dashed line show the curves 
obtained from the combined {J /'4' and ip' data) fit for the temperature depending on which 
parton distribution has been used. 

non-polarization of the final J /ip., there is no correlation between the proton spin and the 
final J /ip spin orientation. All single spin asymmetries are then zero as to the order of 
accuracy of the approximations applied here. 

• The double spin asymmetries for the directly produced J /ip and ip' depend up to a factor 
only on the ratio of the polarized and unpolarized gluon distributions. In case of the 
inclusive J/ip cross section the different mass-scales make the situation more complicated. 

These findings mean a big simplification for the polarized physics because they state that 
the extraction of the polarized gluon density from the double spin asymmetry will not be 
complicated by initial and final state spin correlations. With the model set up in the previous 
chapter we are able to compute the error bars for the double spin asymmetries for J /ip, ip' and 
Xcj production. In the following we collect the expressions for the double longitudinal spin 
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Figure 7: Contribution of various subprocesses to inclusive J /ip production in pA and vrA 
collisions. On the right the ratio of the directly produced J/ip to all J /ip measured is displayed. 
In the middle the ratio between the ip' cross section and the cross section of the directly 
produced J /ip is shown. The shape of this ratio changes with s only due to the evolution of the 
gluon parton density. On the right the ratio XcilXc2 as discussed in the text is displayed. The 
underlying gluon distribution for all figures has been taken from GRV and GRS respectively. 



cross section Actll: 

Actll = — I dxidx2'^T,LL{Sz) (general formula) 

AcTLL(J/^)dir = ^"^^-J dx^ j dx2J's{M]/^)^G{x^,M]/^)^G{x2.Ml 

AaLLix'cj) = i2J+l)^^ldx,ldx2Tp{Ml^,s)^Gix,,Ml^,s)^Gix2,Ml 

AaLL{J/^)inci = A(Tll( J/^)dir + Bt{^' ^ J/^)AaLL{^') + ^ Bt{xcJ ^ J/^)AaLL{x 

J 






XcjeflfJ 



cj) 



(35) 



The - sign in the general formula takes into account that the standard convention for the nu- 
merator of the asymmetry is always anti-parallel spin alignment minus parallel spin alignment. 
For the color singlet (CSM) contributions we find: 



XcO 






^^ \ dx,dx,AG{x,^Ml:)AG{x,^Ml:)6 [l - -^j 



(36) 
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Figure 8: Double spin asymmetry All for inclusive J/ip (left) and ip' (right) hadroproduction. 
For the plot we use the unpolarized parton distribution set GRV and the polarized GSA. For 
the grey error band the assumed luminosity is 0.25 pb~ using a beam polarization of 100%. 

The corresponding unpolarized cross sections can be straightforwardly obtained by replacing 
the polarized by the unpolarized gluon distributions and to remove all minus signs. Then the 
double spin asymmetry All and its statistical error 6 All are simply given by: 



All = SaLL/cr ; 



6A 



Gsa_ 



LL 



Cia, 



O- 



cr± = cr ± Actll . 



(37) 



Fig. 1^ shows the double spin asymmetry for S-wave charmonium production, i.e inclusive J/ip 
and ip' production. For the polarized gluon distribution amplitude we use the leading order set 
gluon A 1^, which we will abbreviate in the following by GSA. For the error band we have 
assumed a luminosity of £ = 0.25 pb~^. It is seen that for inclusive J /ip and ip' production the 
asymmetry at RHIC energies is sizable. It is larger for ip' production, but here also the error 
bars are larger. For P-wave charmonium production (see Fig. P) the asymmetry for the Xc2 
production is partially negative due to the GSM contribution. In case the GSM contribution 
is smaller, i.e. that a smaller value for |i?'^| is more realistic, the asymmetry will go to more 
positive values, so here we see a very fine test of the interplay between GSM and comovers. In 
all cases we get a substantial asymmetry for RHIG energies with small error bars. In general it 
is noticed that the asymmetry decreases with increasing beam energy y/s. This means that in 
addition to the RHIG spin program an polarized experiment like HERA-iV, which is supposed 
to run at ^/s = 40 GeV will also contribute very valuable information for the mechanism how 
charmonium production happens in hadroproduction. One should notice that the asymmetries 
for inclusive J/tp and the Xci production do not depend on the temperature while the other 
asymmetries do, so to measure the asymmetries will give essential new information upon the 
validity of the theory. 

Finally, it has to be stated that the whole calculation still is based on the velocity expan- 
sion p/m within the NRQGD formalism, which is truncated already at the first order. The 
corrections to this may be quite considerable. In this direction we find also the big expansion 
parameter p already determined in f^ and which has been confirmed here. Unfortunately, the 
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Figure 9: Double spin asymmetry All for inclusive Xcj, J = 0,1,2 hadropro duct ion. For the 
plot we use the unpolarized parton distribution set GRV and the polarized GSA. For the grey 
error band the assumed luminosity is 0.25 pb~ using a beam polarization of 100%. 

inclusion of those velocity corrections will destroy the simple relations derived in |^, ^]. It would 
be quite advisable for future studies to develop a formalism, that could test the results of the 
hard comover rescattering from a different stand point which is not based on the approximations 
of the NRQCD. 

7 Summary 

In this work we have tried to describe the measured unpolarized cross section for charmonium 
production through the framework of hard comover rescattering and made predictions for the 
asymmetries in polarized pp scattering. The generic advantage of the hard comover rescattering 
mechanism is that it can explain the non-polarization and the comparatively large value for 
the xi/X2-fraction observed in experiment in a simple and natural way. 

In order to get quantitative results we have expressed the comovers as a thermal cloud of 
gluons. The measured data suggest that about 0.5-1% of the total energy in the collision is 
invested in the formation of the gluon cloud. The fact, that the final state J/ip is unpolarized 
leads us to the conclusion that the displacement parameter b should be small and consistent to 
zero. If this is true then it means that there is no correlation in polarized experiments between 
the initial polarization of the protons and the final polarization of the J/ip and furthermore 
the single-spin asymmetry Sql should be zero, a notion which should be tested by experiment. 
The asymmetries AllH^') and All{Xci) depend only on the ratio of the polarized and unpolar- 
ized gluon distribution amplitudes, while ALL{J/4')mci, ^ll(Xco) and All{Xc2) are in principle 
sensitive to the parameters describing the gluon cloud. The hard comover rescattering picture 
provides an understanding of the formation of onium states in hadroproduction which may give 
the answers to some problems left unsolved in the standard COM and CSM mechanism. The 
RHIC-spin experiment and a possible HERA-A^ will provide very crucial new information upon 
the validity and the consequences of the theory presented here. 
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A Derivation of the xo ^^^ X2 ^(<^?) CSM cross sections 

In this appendix we reproduce the Born cross section for xo and X2 production ((9(a^)). This 
gives a cross check for the formulas derived in [Q and also a cross check for the phase-space 
and flux factors we used in the text to obtain the cross section formulas from the amplitudes. 



The momenta for the gluon fusion amplitude are: 

9i 



■- (m, 0, 0, m) 
Pi = {m, p) 



92 
P2 



(m, 0, 0, — m) 
(m, -p) 



(3^ 
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-ig^v{p2,X) 
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.ni^,-^.+-.V,T^(^ + ^) 



+^,T''{P,-^2 + 
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Pz 



m m^ 



u(,Pi, A) 



+ 7^fahc ei ■ 62 v{p2, A)7^m(pi, A) 
2m 



(39) 



Using now the equation of motion one can write: 
^^T\p^ -i, + mgy^T'' = -T'T^ y^e2 ■ p + i^^si ■ p + m (^i • £273 + i[ei x esj^ToTs) 



(40) 



The spinor combinations can be expressed as follows: 



2mS^^2\ 
2m5^^ 



v{p2,>^h^u{pi,X) 
^(P2,A)7%5m(pi,A) 

v{p2,Xyu{puX) = 2m6l{-2X){e*{2X)-6)V2 



(41) 



Then, using [T"',T'^] = ifabcT'^ and {T°',T^} = dabcT'^, we arrive at the following expression in 
the first order p/m: 



$ 



igUi6)^^[ei x 82] 



1 



-rr^abSij + dabcT^j fabcT- 



iPz 



-2XSlV2iie*i2X) ■ e,){62 ■ p) + (e*(2A) ■ 62)ie, ■ p)) - l—5abS^J + d^bc^ 



m \N, 



Pz ( '^ 
'm\Nr 



+2X5^\i ■ £23 ( —5ab5ij + dabcT,^ 



(42) 
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Hereby we reproduce up to convention dependent phase factors Eq. (4) in 0. For the wave 
function one uses the following expression: 

V'^f ^ (q) = ,Pll. l=e{S,)x--x^Xx ■ (43) 

Now the amplitude for XcJi J = 0, 1, 2 production is given by: 



M^-ep.„j.)^Y.(l ij l)/(w<lw*'^'i"«' 



(44) 



With this we find for the Xco nieson: 



A^(»)(3Po, 0) = -9"j^^^^A^^^,R'l^M' , (45) 

with R'-^ being the first derivative of the quarkonium wave function at the origin, as defined by: 



/f^^-.w-Vii«wi.)- («) 



Then the partonic cross section is given by: 



2^,21 D 12 



The amplitude for the xa meson vanishes identically. For the Xc2 meson we find then: 



M^'^K'P,, J.) = -9'l-ya,^j^,R'Al6]l,, , (48) 

and we obtain for the partonic cross section henceforward: 

^AiAal ^2, ■Jz)cSM - -^^ ^a/j,/2 ' (49) 

B On the alternative of a virtual gluon field T 

A virtual gluon field can be parameterized by the transversality condition: 

r.r;, = g,,, - ^ . (50) 

Taking now the incoming and outgoing charm quark to be on-shell one obtains up to velocity 
corrections: 

ml = {l+py = ml + l'^ + 2p-l^f + 2/om + m^ , (51) 

which results in: 

|roP = l + ^. (52) 
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If IFqI is going to be zero it requires /q = —2m which means that the charm quark gets a 
negative energy, which is unphysical. Now we can investigate how much the incoming charm 
quark needs to be off-shell so that we can work with a real gluon field instead: 

21 

mf. = (l+p) = (l- e)ml + 2p ■ I ^ e ^ — . (53) 

m 

Then, for small enough Iq, we can again treat the off-shellness as a velocity correction along 
the many others we have neglected in the calculation. 
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